Dynamical paths and universality in continuous variables noisy channels 
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We introduce an overcomplete parameter space for two-mode symmetric Gaussian states suitable to address 
the decoherence effects of both Markovian and non-Markovian Gaussian maps. We observe universality of the 
dynamical paths, which do depend only on the initial state and on the effective temperature of the environment, 
with the non Markovianity that manifests itself in the velocity of running over a given path. Universality is also 
seen in the value of discord at the separability threshold and may be exploited to build constants of motions 
valid for both Markovian and non-Markovian maps. We also found that the geometrical constraints provided by 
the structure of the parameter space imply the existence of excluded regions, i.e. sets of states which cannot be 
linked by any Gaussian dynamical map. 
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As soon as entanglement has been recognized as a resource 
for quantum technology, it has been realized that decoherence 
is the main obstacle to overcome. Decoherence appears when- 
ever a system interacts with its environment, so that its dy- 
namics is no longer unitary, but rather described by a noisy 
quantum map, irreversibly driving the system towards relax- 
ation and loss of quantum coherence |fl~||2]. The main effect 
of the interaction with environment is to set up a time scale 
t m over which the dynamics of the system is effectively de- 
scribed by a coarse grained Markovian process towards equi- 
librium. Conversely, for times shorter than t m , the dynamics 
is more involved and the correlations with and within the en- 
vironment play a major role J2K6). In this regime the deco- 
herence may be less detrimental, or even almost unitary, and 
this is why a great attention has been devoted to the study 
of the corresponding non-Markovian maps. Indeed, there are 
evidences that non-Markovian open quantum systems ITl- flOl 
can be useful for quantum technologies ifTTl fl2l . As a con- 
sequence, much attention is currently devoted to the analysis 
of system-environment coupling in order to characterize, con- 
trol, and possibly reduce decoherence in the most effective 
way lfT3l H4l . e.g. by taking advantage of the back-flow of 
information from the environment. 

As a matter of fact, non-Markovian models are more in- 
volved and only few cases can be solved analytically 12 Q3] 
[16). However, these cases are also of great general validity, 
and they can display all of the properties of interest. This is 
especially true for continuous variable systems [17], where 
a set of quantum oscillators excited in a Gaussian state, and 
interacting with their thermal environment, provides an excel- 
lent model for a large class of physical systems in order to 
study non-Markovianity and the decoherence of quantum cor- 
relations. Motivated by the above considerations, we address 
in details the dynamics of quantum correlations between two 
quantum oscillators prepared in a symmetric Gaussian state 



and interacting with local thermal environments. At first, we 
notice that the set of Gaussian states [ 18 1 do not constitute a 
manifold, nor it is convex, and thus geometrical approaches to 
their dynamics are not considered particularly appealing. At 
variance with this belief, we address the study of decoherence 
by representing dynamical paths in a suitable, overcomplete, 
parameter space, involving entanglement, Gaussian discord 
and the overall purity of the state. These variables proved to 
offer a suitable framework to compare non-Markovian maps 
and their Markovian counterparts, and to show which proper- 
ties do, and notably do not, distinguish Markovian and non- 
Markovian processes. In particular, we observe universality of 
the dynamical paths, which do not depend on the specific fea- 
tures of the environment spectrum and are determined only 
by the initial state and the effective temperature of the envi- 
ronment. The non-Markovianity of the system only changes 
the velocity of running over a given path. This behavior al- 
lows ones to map non-Markovian processes onto Markovian 
ones and it may reduce the number of parameters needed 
to study a dynamical process, e.g. it may be exploited to 
build constants of motions valid for both Markovian and non- 
Markovian maps. Universality is also observed in the value 
of discord at the separability threshold, which moreover is 
a function of the sole initial conditions in the limit of high 
temperature. Finally, we find that the geometrical constraints 
provided by the structure of the parameter space implies the 
existence of excluded regions, i.e. sets of states which cannot 
be linked by any Gaussian dynamical map. 

Dynamical evolution — We address the dynamical decoher- 
ence of two oscillators of frequency uio, each coupled to its 
own bosonic environment made of modes at frequencies uik- 
The baths are separated and of equal structure, such that the 
evolution preserves any initial symmetry between the two os- 
cillators. The system-bath interaction Hamiltonian is given 
by Hi = aY, k j k (X 1 q 1 ^ k + X 2 q 2 ,k) where a is the cou- 
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pling, the complex {jk}k=i,2.. modulate the dispersion over sinh 2 r(2j/ T + 1) + v T , while the c coefficients describe the 



the bath's modes, and Xi — (Sj + gJ)/2, i — 1,2 and 
Qk.i = (bk.i + b\ ■ )/2 denote the canonical operators of the 
systems' and baths' modes respectively. In the weak cou- 
pling limit a <C loq, and discarding the counter rotating terms, 
we get a non-Markovian master equation for the dynami- 
cal evolution of the density operator g describing the quan- 
tum state of the two oscillators. Assuming that the initial 
state is a two-mode Gaussian state with zero amplitude the 
master equation may be transformed into an evolution equa- 
tion for the two-mode covariance matrix lfT9l l20l (CM) a, 
a jk = iTT[g{R j R k +R k R j )],K={X 1 ,Y 1 ,X 2 ,Y 2 }, whose 
solution may be written as 



e- r (*Vo + A r (i) 



(1) 



where cto is the initial CM, 

r(t) = fds^s), A r (t) = e- r ^ f ds e r(s) A(s) 
Jo Jo 

and where, upon defining the spectrum of environment as 
j(ui) = Tj k \j k \ 2 8{uj — ui k ), the diffusion and dissipation func- 
tions are given by Ell : 

r-t />00 

A(t) — a 2 ds dui j(uj) coth(a>/3/2) cos(ws) cos(wos) 
Jo Jo 

j(t) = a 2 / ds duj sin(ws) sin(wos) (2) 
Jo Jo 

where (3 = 1/ kT. At high temperatures the damping is negli- 
gible and the first coefficient is dominant, while at lower tem- 
peratures they have the same order of magnitude. 

The non-Markovian features are embodied in the time de- 
pendence of the coefficients Ap(i) and T(t), describing dif- 
fusion and damping respectively, which for times t < r M are 
strongly influenced by the whole spectrum of the environment 
ll22l |23l . On the other hand, for times t ^> t m the coeffi- 
cients achieve their Markovian limiting values. In particu- 
lar we have lim.{_>. +00 7(f) = a 2 |j(oj )| 2 = j M , such that 
I\t) = j M t, A r (t) = (1 - e'" lMt ){2n T + 1), and the solu- 
tion ([!]) rewrites as a(t) = e~ 7Mt cr + (1 - e _7M *)er T , where 
(7 T = [n T + i)I is the CM of the stationary state, i.e. a state 
in thermal equilibrium with the environment and a population 
of n T = {er u — thermal photons. 

Symmetric Gaussian states — Symmetric Gaussian states 
are those with a CM that can be recasted (via local opera- 
tions) in a form depending on two real parameters a(a, c) = 
al4+ccri(8)(73, the er^ 's being Pauli matrices. This correspond 
to preparing the two oscillators in a squeezed thermal state, 
with density operator of the form g(r, )n T = S(r)v®vS* '(r), 
where v is a single mode thermal state with v T photons and 
S(r) = e r ( a i a 2+aia 2 ) is the two-mode squeezing opera- 
tor. The link with the parameter of the CM is given by 
a = (y T + i)cosh(2r) and c = (v T + i)sinh(27-). The 
diagonal elements a = | + n depends on the (equal) pop- 
ulation (mean photon number) of the two subsystems n = 



correlations among them. Uncertainty relations impose a con- 
straint which reads E4ll \a — c| > \. 

As a matter of fact, the representation in terms of the coef- 
ficients a and c does not fully illustrate the correlation prop- 
erties of a state. In particular, it does not allow to analyze 
the relations between different kinds of quantum correlations, 
as entanglement or discord, in a dynamical context. To this 
aim we introduce a different (overcomplete) parametrization 
involving the overall purity of the state /i = Tr[g(r, n T ) 2 ], its 
entanglement expressed in terms of the minimum symplectic 
eigenvalue A = a — c (the state is separable iff A > ^) and 
the Gaussian quantum discord, which for symmetric Gaussian 
states, may be written as Il25ll26l 

_ 2e 2 
1 + 2a 



D{a,c) = h(a) - 2h 



(3) 



where h(x) = (x + ^)log(a; + ^) — (x — ^)log(x — ^). 
Of course, the resulting parameter space is overcomplete and 
the third parameter is a function of the other two E7\ . In 
the following, we will describe the dynamical paths in the 
three-dimensional space (/i, A, D). These will lie over the sur- 
face individuated by the constraint D = D(/i, A), and in the 
region satisfying the uncertainty relations, which rewrite as 
^<|A- 2 . 

Markovian dynamics — The Markovian master equation de- 
pends on the environment (effective) temperature and on the 
damping j M , nonetheless the Markovian dynamical paths de- 
pend exclusively on the (effective) temperature of the envi- 
ronment. The damping only affects the speed of running over 
a dynamical path, but not its shape, and the rate c(t) /co = 
e -7M* determines in a unique way the rate a{t)/a^. In the 
left panel of Fig. [T]we report Markovian paths for different 
temperatures, assuming the two oscillators initially prepared 
in a two-mode squeezed vacuum (TWB) g(ro, 0), i.e a pure 
maximally entangled state. As it is apparent from the plot, 
two limiting paths emerge at low and high temperatures. The 
transition from one regime to the other occurs continuously 
by raising the temperature, and we see that the high temper- 
atures limit is already achieved for temperatures leading to 
n T j sinh 2 (ro) > 3. Two other phenomena are revealed by 
this representation: (i) the value of the discord at the separa- 
bility threshold (A = |) depends only on the initial squeezing 
ro and approaches a universal curve in the high temperature 
limit; (ii) for a given initial state g(ro,0) there are STS that 
cannot be reached during any Markovian decoherence pro- 
cess, despite the fact that they have reduced entanglement and 
purity compared to the initial state. 

Non Markovian dynamics — As mentioned in the introduc- 
tion, non-Markovian dynamics may display remarkable dif- 
ferences from their Markovian counterpart during the initial 
transient when t < t m . Entanglement oscillations may oc- 
cur, and the separability threshold may be delayed or acceler- 
ated depending on the spectrum of the environment. A ques- 
tion thus arises on whether these differences also affect sig- 
nificantly the dynamical path in the space of parameters. As 
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FIG. 1: Left: Markovian thermalization paths at different tempera- 
tures. The initial state is a STS with ro = 1.2, the different paths 
correspond to n T = 0, 0.1, 0.5, 1, 10. Solid blue line correspond to 
thermal states with zero discord, solid red line denotes the separa- 
bility threshold. The high-temperature limit is already achieved for 
n T > 3. Right: dynamical paths in the high-temperatures limit for 
different values of initial entanglement. The curves correspond to 
r = 0.5,0.7, 1,1.5,2. 

we will see, the answer is negative, and universality occurs. 
The results about the dynamics that we are going to discuss 
are independent of the particular choice of environment spec- 
trum. However, in order to show some numerical solutions, 
we employ few examples corresponding to white noise and 
to both Ohmic and super-Ohmic spectral densities with cut- 
off lu c [28 1 . Let us start by analyzing the high temperature 
regime, where over a timescale r ~ t m we can neglect the 
damping T(t) (it becomes relevant over times r ~ j~ ^> t m , 
which is definitely in the Markovian regime). Short time non- 
Markovian dynamics is thus due to the behavior of the heating 
function Ar (t) and, in turn, is very sensitive to the details of 
the environment spectrum j(co). In this limit non-Markovian 
effects can be seen during the whole decoherence process, 
with entanglement oscillation across the separability threshold 
|[T9l . The dynamics is driven by the approximate dynamical 
equation 

a t ~a Q + J ds AO) -| (4) 

corresponding to a(t) = ao + \ J Q dr A(t) and c(t) = cq. 
The last condition imposes a constraint to the dynamical 
paths, which is the same independently on whether the dy- 
namic of a(t) is Markovian or displays oscillations, as long as 
a(t) > ao Vf and a(t) — » a T . In other words, the paths are the 
same of the Markovian case, and the possible oscillations of 
a(t) only influences the speed of running over the dynamical 
path. In the right panel of Fig.[T]we show the dynamical paths 
for different values of the initial squeezing r . 

The condition c(i) = cq also implies that the Gaussian 
discord may be written as D(a,c) ~ D(X + cq,cq), i.e. 



it depends on the temperature and on the initial squeezing 
ll29l . At the separability threshold we have D sep = -D(i[l + 
sinh(2ro)], ro) = /( r o), i.e. the discord at separability is a 
universal function of the initial squeezing. In Fig|2]we show 
the Gaussian discord at separability as a function the initial 
squeezing. The solid black line correspond to the high temper- 
ature approximation /(ro), whereas the colored dots are ob- 
tained using the full non-Markovian solutions, obtained taking 
into account the damping and different environment spectra. 
We first point out that there is excellent agreement among the 
two solutions already for temperature about n T < 1 and inde- 
pendently on the environment spectrum. We also notice that 
D sep saturates to a limiting value d* = lim^-j-oo + 
sinh(2r )],r ) = -1 + 2 log 2 ~ 0.3863 when the initial 
squeezing increases. The initial squeezing needed for achieve 
the saturation regime increases with temperature. As it may 
be seen from the plot, for high temperatures, i.e. for n T > 1, 
it is about ro ~ 2. 




FIG. 2: Discord at separability threshold as a function of the ini- 
tial squeezing. Solid black line represents the universal function 
/(ro), whereas the dots are the numerical solutions of the full non- 
Markovian dynamics in Eq. |TJ for n T = 10 and for three differ- 
ent environment spectra corresponding to ohmic, superohmic, and 
white noise spectral density respectively (red circles, green dia- 
monds, and blue squares). The dashed line is the high-temperature 
high-squeezing limiting value d* ~ 0.3863. The gray lines de- 
note the Markovian curves D(| + c(t Bep ), c(t sep )) for n T = 
0.5, 0.1, 10 -2 , 10 -3 respectively. We also report the numerical solu- 
tions of the full non-Markovian dynamics for n T = 0.5, 0.1 and the 
same three spectra, whereas for n T = 10~ 2 , 10~ 3 the separability 
threshold is definitely in the Markovian regime. 

The dashed gray lines in Fig. |2]denote the Gaussian discord 
at separability D sep — + c(t sep ), c(t sep )) for various 
(low) temperatures, whereas the dots are obtained using the 
full non-Markovian solutions for different environment spec- 
tra. At low temperatures the damping j(t) and the heating 
function A(t) become of the same order of magnitude and 
thus the separability threshold t sep does depend on the envi- 
ronment spectrum. On the other hand, separability is always 
achieved in the Markovian regime, and thus D sep is a univer- 
sal quantity. For times t < t m , there is a competition be- 
tween 7(f) and A(t) and in principle, one would not expect a 
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universal behavior. However, low temperature and weak cou- 
pling make the effect of damping and heating very weak, with 
appreciable perturbation of the initial state only after a long 
time. In other words, any dynamical effect of the interaction 
is taking place in the Markovian regime, thus being univer- 
sal and independent on the environment spectrum. This also 
means that the dynamical paths in the left panel of Fig. [T] le- 
gitimately describe non-Markovian dynamical trajectories at 
low temperatures. More generally, any path-dependent prop- 
erty can be checked analytically using the set of Markovian 
equations and then extended to the non-Markovian regime, 
where an analytic approach would be unfeasible. In particu- 
lar, let us introduce the rescaled time r = Ft, and recall that 
in the Markovian regime we have d T X = e~ r (A T — Ao) and 
d T (Xfi)^ 1 = e~ T [(Ao/io) -1 + 4A T ], where pedices 0/T re- 
fer to initial/stationary state. Then, any constant of motion, 
e.g. C = A + y/(4A/x), with y = (A T — A )/A r + (4^ Ao) _1 
built using the Markovian dynamical equation is a constant 
of motion also in the non-Markovian regime, independently 
on the environment spectrum. The temperature dependence 
disappears in the high temperatures limit. 

Discussion and conclusions — We have introduced an over- 
complete parameter space for two-mode symmetric Gaus- 
sian states suitable to address the decoherence effects of both 
Markovian and non-Markovian Gaussian maps, and to reveal 
universality of the dynamical paths. We have shown that 
Markovian paths coincide with non-Markovian ones and thus 
all path-dependent properties are not depending on the en- 
vironment spectrum, e.g. it is possible to build constant of 
motions valid for both Markovian and non-Markovian maps. 
Dynamical paths do depend only on the initial state and on the 
effective temperature of the environment, and the non Marko- 
vianity behavior is only seen he velocity of running over a 
given path. Our results also show that a non-Markovianity 
measure ll7l [T0l l cannot be based on any geometrical charac- 
terization of the dynamical paths. 

Universality is also observed for the value of discord at the 
separability threshold, which moreover depends on the sole 
initial squeezing in the high temperature limit. We also found 
that the geometrical constraints provided by the structure of 
the parameter space implies the existence of excluded regions, 
i.e. sets of Gaussian states which cannot be linked by any 
Gaussian dynamical map. Finally, we emphasize that univer- 
sality of dynamical paths is related only to the assumptions 
of weak coupling and linear interaction and therefore it may 
be conjectured that it represents a more general feature, char- 
acterizing any open quantum system admitting a Markovian 
limit. 
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